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(ABSTRACT) 

^^This  dissertation  considers  the  discrete-time  control  of 
a  spacecraft  consisting  of  a  rigid-^platform  with  retargetable 
flexible  antennas.  The  mission  consists  of  independent 
minimum-time  maneuvers  of  each  antenna  to  coincide  with  pre¬ 
determined  lines  of  sight,  while  the  platform  is  stabilized 
in  an  inertial  space  and  elastic  vibration  of  the  antennas  is 
suppressed.  The  system  is  governed  by  a  set  of  linearized, 
time-varying  equations  of  motion.  A  discrete-time  approach 
permits  consideration  of  the  time-varying  nature  of  the  system 
in  designing  the  control  law. 

Both  global  and  decentralized  controls  are  proposed  for 
a  noise-free  system  with  full-state  feedback.  Initially,  a 
time-varying  linear-quadratic  regulator  (LQR)  is  implemented, 
followed  by  two  types  of  decentralized  controllers.)^  First, 
a  collocated  control  law  is  devised  in  which  actuator  forces 
are  based  on  the  position  and  velocity  at  the  actuator  locations. 
Next,  a  new  method  called  Substructure-Decentralized  Control 


is  proposed,  where  each  flexible  substructure  is  controlled 


based  on  state  measurements  associated  with  the  substructure 
modes  of  the  separately  modeled  appendages. 

In  both  global  and  decentralized  cases,  a  linear  control 
law  is  first  implemented  coupled  with  an  open-loop 
disturbance-accommodating  control  based  on  the  known  inertial 
disturbances  caused  by  the  maneuver.  Elastic  motion  is  next 
controlled  using  nonlinear  (on-off)  antenna  controllers  for 
each  decentralized  case.  For  Substructure-Decentralized 
Control,  the  controls  translate  into  quantized  actual  controls. 
Lastly,  nonlinear  (on-off)  control  laws  are  also  used  to  control 
the  rigid-body  motion  for  each  case. 

Next,  the  problem  of  controlling  the  time-varying  system 
in  the  presence  of  noisy  actuators  and  sensors  is  examined. 
It  is  assumed  that  only  displacements,  not  velocities,  are 
sensed  for  both  rigid-body  and  elastic  motion,  making  state 
reconstruction  also  necessary.  A  discrete-time,  full-order 
Kalman  filter  is  constructed  for  the  time-varying  system.  A 
pseudo-decentralized  control  is  proposed  whereby  feedback 
controls  are  based  on  system  state  estimates.  As  before,  both 
linear  and  nonlinear  controls  are  implemented.  For  each  case 
mentioned,  a  numerical  example  is  presented  involving  a 
spacecraft  with  a  single  flexible  maneuvering  antenna. 
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1  INTRODUCTION 


As  the  complexity  of  future  space  systems  and  missions 
increases,  e.g.,  NASA's  Space  Station,  the  problem  of  maneuver 
and  control  of  flexible  spacecraft  takes  on  added  importance. 
The  orientation  of  simple  flexible  spacecraft  along  a  desired 
line  of  sight  has  been  generally  addressed  by  slewing  the 
entire  vehicle  and  suppressing  vibration  in  the  flexible 
components  during  and/or  after  maneuver  completion.  However, 
instances  arise  in  which  it  is  more  efficient  to  reorient 
specific  components,  e.g.,  an  antenna,  relative  to  a  main  body 
while  stabilizing  the  main  body  in  an  inertial  space  and 
suppressing  elastic  motion  of  the  flexible  parts  of  the 
spacecraft.  Such  an  approach  becomes  particularly  attractive 
in  the  case  of  spacecraft  consisting  of  a  main  rigid-body  and 
several  flexible  substructures,  each  requiring  independent, 
simultaneous  retargeting.  This  dissertation  addresses  the 
control  of  such  a  system. 

1 . 1  Background 

Since  the  beginning  of  the  "Space  Age"  more  than  30  years 
ago,  the  dynamics  and  control  of  flexible  spacecraft  have  been 
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a  growing  concern  (Ref.  39)  .  The  problem  of  efficient  derivation 
of  the  equations  of  motion  for  complex  spacecraft  has  been  a 
constant  goal  in  the  analysis  of  systems  such  as  the  one 
considered  in  this  dissertation.  Indeed,  equations  of  motion 
for  flexible  spacecraft  are  typically  very  complicated,  leading 
to  many  recent  attempts  to  develop  more  convenient  expressions 
for  spacecraft  motion.  Several  of  these  are  compared  in  Ref. 
16.  Another  effective  modeling  technique  using  Lagrange's 
equations  of  motion  for  flexible  bodies  in  terms  of  quasi- 
coordinates  (Ref.  23)  was  proposed  by  Meirovitch  (Ref.  33) . 

As  the  complexity  of  missions  has  grown  beyond  simple 
station  keeping  to  include  maneuvering,  control  design  has 
also  become  more  difficult.  The  control  of  large  flexible 
structures  in  general  is  treated  in  several  works,  e.g.  Refs. 
2,4,15,16,22,25-27,51.  Studies  involving  single-axis  maneu¬ 
vers  of  simple  spacecraft  consisting  of  a  main  rigid-body  and 
a  number  of  flexible  appendages  with  fixed  orientation  relative 
to  the  main  body  provided  important  insights  into  problems 
confronting  engineers  (Refs.  5,7,8,30,31,42,44,47,48).  Turner 
and  Junkins  (Ref.  47) ,  Turner  and  Chun  (Ref.  48)  and  Breakwell 
(Ref.  8)  proposed  methods  for  simultaneous  maneuver  and 
vibration  suppression  of  such  systems,  while  Baruh  and  Sil- 
verberg  (Refs.  5  and  7)  separated  the  maneuver  and  vibration 
suppression  tasks.  More  recent  work  has  centered  on  NASA's 
Spacecraft  Control  Laboratory  Experiment  (SCOLE)  and  the 
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possibility  of  missions  involving  control  of  flexible  struc¬ 
tures  carried  into  orbit  by  the  space  shuttle  (Refs.  29,42,49)  . 
The  objective  of  SCOLE  is  to  reorient  the  structure  line  of 
sight  in  minimum  time  with  limited  control  authority.  As  in 
the  case  of  simpler  systems,  the  SCOLE  is  concerned  with 
reorienting  the  entire  spacecraft. 

For  more  complex  systems  involving  rotational  maneuvering, 
the  derivation  of  the  equations  of  motion  by  means  of  Lagrange's 
equations  for  flexible  bodies  in  terms  of  quasi-coordinates 
is  particularly  attractive.  The  beauty  of  this  method  is  that 
the  equations  of  motion  are  based  on  body-fixed  coordinates, 
hence  it  is  convenient  to  design  a  control  law  based  on  these 
body-fixed  coordinates.  This  approach  was  successfully  applied 
by  Meirovitch  and  Kwak  (Ref.  35)  to  multi-body  systems  in  which 
individual  substructures  are  reoriented  to  coincide  with 
desired  lines  of  sight.  Figure  1  shows  a  typical  spacecraft 
of  the  type  described  above  consisting  of  a  rigid  main  body 
and  several  flexible  appendages. 

The  control  design  problem  for  a  system  of  the  type  shown 
in  Fig.  1  was  first  treated  by  a  continuous-time  approach  by 
Meirovitch  and  Kwak  (Ref.  32)  under  the  assumption  that  the 
time-varying  terms  in  the  coefficients  are  sufficiently  small 
that  they  can  be  ignored  for  the  purpose  of  control  design. 
This  is  based  on  the  premise  that  the  appendages  are  relatively 
small  and  that  the  maneuvers  are  "slow."  Reference  32  also 
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assumed  noise-free,  continuous,  full-state  feedback  by  means 
of  perfect  actuators  and  that  each  antenna  maneuver  is  carried 
out  open-loop  using  a  bang-bang  control  law,  while  platform 
inertial  stabilization  and  appendage  vibration  suppression  are 
performed  simultaneously.  Reference  19  contains  a  summary  of 
Refs.  32,  34,  and  35  that  propose  a  disturbance-minimization 
technique,  a  proportional-plus-integral  (PI)  control  and  a 
perturbation  method,  respectively,  to  control  the  spacecraft. 
A  perturbation  method  was  applied  to  the  time-varying  system 
in  Ref.  34  using  both  an  "adiabatic  approximation"  (Ref.  12) 
and  an  integral  method  to  update  the  system  parameters  and 
calculate  feedback  control  gains  at  discrete  intervals  during 
a  maneuver  while  employing  a  continuous  time  controller.  Gains 
calculated  using  both  techniques  are  held  constant  between 
updates . 

The  time-varying  nature  of  the  system  was  considered  by 
Meirovitch  and  France  (Ref.  36)  using  a  discrete-time  approach. 
The  resulting  control  laws  assumed  noise-free,  full-state 
feedback  and  perfect  actuators.  Both  decentralized  and  global 
linear  control  laws  were  demonstrated,  as  well  as  a  combination 
of  linear  and  nonlinear  decentralized  controls  to  stabilize 
the  rigid-body  motions  and  elastic  deformations,  respectively. 
Reference  38  extended  the  work  of  Ref.  36  by  considering  noisy 
actuators  and  sensors,  where  the  sensors  measure  elastic 
deformations  of  the  antennas  relative  to  local  coordinate 
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frames.  This  sensor  output  is  then  processed  by  a  Kalman 
filter  that  provides  current  estimates  of  the  system  states, 
upon  which  each  antenna's  pseudo-decentralized  control  input 
is  based. 

The  use  of  strain  gauges  in  measuring  elastic  motions  has 
been  discussed  in  several  papers  (Refs.  9,15,18,41)  .  Reference 
18  demonstrated  that  strain  gauges  could  be  used  to  measure 
very  high  frequency  phenomena.  Optical  sensors  also  show  great 
promise  in  low-noise  sensing  of  elastic  motions  (Refs.  43,49) 
as  well  as  in  distributed  sensing  (Refs.  40,50). 

Nonlinear  control  of  flexible  structures  using  quantized 
controls  was  proposed  in  Ref.  28,  while  on-off  attitude  control 
of  spacecraft,  and  rigid-body  systems  in  general,  is  discussed 
in  many  sources  (e.g..  Refs.  1,17,46).  Several  decentralized 
control  techniques  for  control  of  large  flexible  structures 
have  also  been  proposed  (Refs.  3,22,45).  Uniform  Damping, 
developed  by  Silverberg  (Ref.  45)  is  attractive  because  control 
forces  at  discrete  actuator  points  can  be  determined  inde¬ 
pendently  of  stiffness  properties  given  position  and  velocity 
information  at  the  actuator  location.  A  method  of  substructure 
control  based  on  the  Component  Mode  Synthesis  method  has  been 
proposed  by  Young  (Ref.  51)  .  Reference  36  demonstrates  the 
use  of  Uniform  Damping  in  a  discrete-time  formulation  for  the 
system  in  question,  and  it  proposes  a  Substructure- 
Decentralized  Control. 
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The  design  of  both  continuous-time  and  discrete-time 
full-order  Kalman  filters  to  estimate  states  in  a  system  with 
incomplete  and  noisy  measurements  is  well  known  (Refs. 
11,13,20,27).  In  the  case  in  which  some  measurements  are 
noise-free  and  others  are  not,  i.e.,  a  singular  sensor  noise 
covariance  matrix,  some  attempts  have  been  made  to  define  a 
reduced  order  Kalman  filter  (Refs.  10,21),  although  the  com¬ 
putational  savings  when  applied  to  real  systems  remains  in 
question. 

Much  of  this  dissertation  is  based  on  the  work  presented 
in  Refs.  36  and  38.  It  extends  this  earlier  work  by  including 
the  use  of  a  nonlinear  (bang-off -bang)  control  law  to  stabilize 
the  rigid-body  motions  and  by  considering  a  pseudo-col located 
control  to  suppress  vibration  of  the  flexible  appendages.  As 
before,  each  additional  approach  is  demonstrated  by  means  of 
a  numerical  example  involving  a  spacecraft  with  a  single 
flexible  antenna  undergoing  a  45°  reorientation  relative  to 
the  platform,  while  the  platform  is  stabilized  relative  to  an 
inertial  space.  Results  for  the  various  cases  are  then  compared 
and  discussed. 

1.2  Organization  of  tha  Dissertation 

Chapter  2  gives  a  summary  of  the  equations  of  motion  for 
a  spacecraft  consisting  of  a  main  rigid-body  and  several 
independently  retargetable  flexible  antennas  developed  in  Ref. 
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26.  Coefficient  matrices  for  these  equations  are  shown  in 
detail  in  Appendix  B. 

Chapter  3  presents  both  global  and  decentralized  control 
laws  for  the  deterministic  (noise-free)  system  assuming 
full-state  feedback.  First,  a  discrete-time  version  of  the 
linear  disturbance-accommodating  control  derived  in  Ref.  32 
is  presented  based  on  the  inertial  disturbances  caused  by  the 
known  antenna  maneuver.  Then,  a  step-varying  linear  quadratic 
regulator  is  proposed.  Both  linear  and  nonlinear  collocated 
controllers  are  discussed,  with  the  linear  case  based  on  the 
concept  of  Uniform  Damping.  A  new  approach,  Substructure- 
Decentralized  Control,  is  then  presented  using  both  linear  and 
nonlinear  control. 

The  problems  associated  with  controlling  a  noisy  system 
with  incomplete  measurements  is  addressed  in  Chapter  4.  First, 
the  unique  problems  associated  with  sensor  selection  and  output 
processing  are  discussed,  as  are  the  advantages  and  disad¬ 
vantages  of  several  specific  sensor  types.  Design  of  an 
appropriate  step-varying,  full -order  Kalman  filter  is  then 
presented,  followed  by  a  discussion  of  the  effects  of  such  a 
global  estimation  on  controller  design.  Next,  the  concept  of 
pseudo-decentralized  control  is  addressed  for  each  of  the 
decentralized  control  techniques  presented  in  Chapter  3. 

Chapter  5  presents  numerical  examples  of  each  control  law 
considered  for  both  the  deterministic  and  stochastic  systems. 


8 


Each  case  is  evaluated  on  a  system  consisting  of  a  main 
rigid-body  and  one  flexible  antenna  undergoing  a  45° 
reorientation  relative  to  the  platform,  while  the  platform  is 
stabilized  relative  to  an  inertial  space.  Finally,  in  Chapter 
6,  the  work  contained  in  this  dissertation  is  summarized  and 
the  results  of  the  numerical  examples  are  compared  and  dis¬ 


cussed. 


2  EQUATIONS  OF  MOTION 


The  equations  of  motion  for  the  system  shown  in  Fig.  1 
were  derived  in  Ref.  32  using  a  new  formulation  of  Lagrange's 
equations  for  flexible  bodies  in  terms  of  quasi-coordinates 
(Ref.  33) .  Following  is  a  brief  summary  of  these  equations. 
Referring  to  Fig.  1,  we  identify  a  set  of  inertial  axes  XYZ  , 
a  set  of  body  axes  xyz  attached  to  the  rigid  platform  and  a 
set  of  body  axes  x.y.z,  embedded  in  the  typical  flexible 

appendage  e,  (e=  1,2 . N)  .  The  position  vector  of  a  point  in 

the  rigid  body  is  given  by  £r  =  £0  +  r  and  that  in  the  appendage 

e  can  be  written  as  Ra  -  R0  +  ro<  +  r#  +  ua ,  (e=l,2 . A/),  where  Ro 

is  the  radius  vector  from  0  to  o,  r  is  the  position  vector  of 
a  point  in  the  rigid  body  relative  to  xyz  ,  ro<  is  the  radius 
vector  from  o  to  e,  rg  is  the  position  vector  of  a  nominal  point 
in  the  undeformed  appendage  relative  to  x.y.z,  and  ue  is  the 
elastic  displacement  of  that  point.  Vector  R0  is  expressed 
in  terms  of  components  along  XYZ  ,  r  and  ro  in  terms  of  components 
along  xyz  and  r#  and  ua  in  terms  of  components  along  x,y#ze. 
The  velocity  vector  of  o  can  be  written  in  terms  of  components 
along  xyz  in  the  form  V_o  =  C  Ro ,  where  C  is  the  matrix  of  direction 
cosines  between  xyz  and  XYZ  and  Ro  is  the  velocity  vector  of 
o  in  terms  of  components  along  XYZ  .  The  angular  velocity 
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vector  of  axes  xyz  in  terms  of  components  along  xyz  is  given 
by  to-£>0,  where  0  is  a  vector  of  angular  velocities  0j  and  D 
is  a  matrix  depending  on  the  angular  displacements  0<  (i  =  1,2,3). 
In  view  of  the  above,  the  velocity  vector  of  a  point  on  the 
rigid  body  in  terms  of  components  along  xyz  is 

Kr«K/u)Xr  (2-1) 

and  that  of  a  point  in  appendage  e  in  terms  of  components  along 
x.y.z,  is 

Ke  =  Ee(Ko  +  ^x£o<,)  +  (Ee^  +  ^e)x(Ce  +  “B)  +  ^e'  e-1,2 . N  (2-2) 

where  coe  is  the  angular  velocity  of  axes  x„y,z„  E„  is  a  matrix 

of  direction  cosines  between  axes  x.y.z,  and  xyz  and  ve  is  the 
elastic  velocity  of  the  point  in  the  appendage  relative  to 
x,y,z„  ue  =  u#.  For  the  proposed  maneuver,  the  angular  velocity 
vectors  u>e  of  x,yez9  relative  to  xyz  are  given,  so  that  the 
rotational  motions  of  the  appendages  relative  to  the  platform 
do  not  add  degrees  of  freedom  to  the  system. 

The  equations  describing  the  rigid-body  translations  and 
rotations  of  the  system  are  ordinary  differential  equations 
and  those  for  the  elastic  motions  of  the  appendages  are  partial 
differential  equations,  so  that  the  equations  of  motion  are 
hybrid.  Because  control  design  for  systems  described  by  hybrid 
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equations  is  not  feasible,  we  must  discretize  the  system  in 
space.  This  is  accomplished  by  expressing  the  elastic  dis¬ 
placements  as  linear  combinations  of  space-dependent  admissible 
functions  multiplied  by  time-dependent  generalized  coordinates 
of  the  form 


“.(£.•*)  “<t>.(C,)2. (0.  e  =  l  ,2 


.,/V 


(2-3) 


where  4>e  is  a  matrix  of  admissible  functions  and  qe  is  a  vector 

of  generalized  coordinates.  The  equations  of  motion,  dis¬ 
cretized  in  space,  are  given  in  Ref.  32  and  will  not  be  repeated 
here  in  full.  Assuming  small  motions  of  the  stabilized  rigid 
platform,  a  linearized  version  of  these  equations  has  the  state 
form 


x(0  =  A(0x(0  +  B(0/(0  +  D(0d(0 


(2-4) 


where  x(t)  =  [^QTqJq^...qJ/V^a)TqJq^...q^/]T  is  a  state  vector,  in 

which  9  is  a  symbolic  vector  of  angular  displacements  of  the 
platform  and 


A(0« 


_0 _ I _ J_ 

-‘(0K(0  r-M~(7)G(o; 


(2  -  5a) 


0  0 
B(0“  “  -7 - .  -  ‘  D(0-  -  -  , - 

M  ‘(OB  (Q  M  1  (0 


(2-56, c) 
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are  coefficient  matrices.  Moreover,  d(0  is  a  vector  of 

disturbances  caused  by  the  maneuver  and  /(0  =  [F_’J A4*T/[, 

l\nlfT2ilT22--'ir2n2Lli-'-lTNnJT  is  the  control  force  vector,  in  which 
F_*0  and  are  actuator  force  and  torque  vectors  acting  on  the 
rigid  platform  and  /  are  the  actuator  forces  acting  on 
appendage  e  at  point  ,  where  appendage  e  has  na  actuators. 
Explicit  expressions  for  the  matrices  M(f)/  ^(0/  G( 0  and  B(0 
are  given  in  Appendix  A. 


3  DETERMINISTIC  DISCRETE-TIME  CONTROL 


3 . 1  Introduction 

For  the  system  considered  in  this  dissertation,  maneuvering 
of  each  antenna  is  carried  out  open-loop  using  a  bang-bang 
control  law,  while  simultaneously  performing  platform  stabi¬ 
lization  and  appendage  vibration  suppression.  To  this  end, 
the  control  law  design  process  begins  by  considering  a 
discrete-time  approach  that  permits  consideration  of  the 
time-varying  nature  of  the  system.  The  resulting  step-varying 
controller  can  then  be  applied  to  the  time-varying  spacecraft. 
It  is  initially  assumed  that,  as  in  Refs.  32  and  34,  the  sensors 
and  actuators  are  noise  free,  i.e.,  that  the  system  is  det¬ 
erministic,  and  all  the  states  necessary  to  compute  the  feedback 
controls  are  available.  Therefore,  no  state  estimation  or 
filtering  is  required. 

In  this  section,  the  control  process  is  divided  into  open- 
and  closed-loop  components  for  linear  control.  The  open-loop 
control  is  a  the  discrete-time  adaptation  of  the  continuous-time 
disturbance-accommodating  controller  presented  in  Ref.  32. 
Next,  both  global  and  decentralized  control  laws  are  considered. 
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The  global  control  is  a  discrete-time  linear  quadratic  regu¬ 
lator.  Two  different  decentralized  approaches  are  discussed: 
Colocated  Control  and  Substructure-Decentralized  Control.  In 
the  case  of  decentralized  control,  both  linear  and  nonlinear 
controllers  are  proposed  for  rigid-body  stabilization  and 
vibration  suppression. 

3.2  Disturbance-Accommodating  Control 

We  propose  to  control  the  system  by  means  of  a  combination 
of  open-  and  closed-loop  controllers,  taking  into  account  the 
time-varying  nature  of  the  coefficient  matrices  and  the  fact 
that  for  an  open-loop  maneuver  the  inertial  disturbances  are 
known.  Hence,  the  control  is  separated  into  open  and  closed-loop 
components  as  follows: 

/(0-/o(0  +  /c(0  (3-1) 

Then,  with  reference  to  Eq.  (2-4),  the  open-loop  control  is 
chosen  to  satisfy 

B(O/o(O  +  D(Od(O“0  (3-2) 

Considering  Eqs.  (2-5b,c),  then 


B‘(O/o(O  +  d(t)  =  0 


(3-3) 


15 


so  that  the  open-loop  disturbance-accommodating  control  is 

given  by 

/o(0  =  -[B*(0]fd(0  (3-4) 

where 

[B*(0]t  =  [(B*(0)TB*(0]‘,[B*(0]T  (3-5) 

is  the  pseudo-inverse  of  B*(0. 

In  the  case  of  a  discrete-time  system,  the  open-loop 
control  at  step  k  becomes 

/<,(*)  =  -[  B*(fc)]td(A:)  (3-6) 

where  each  of  the  terms  in  Eq.  (3-6)  are  evaluated  at  time 
t  =  kT  as  described  in  the  next  section.  This  open-loop  control 
can  be  implemented  in  conjunction  with  any  linear  closed-loop 
control . 

3 . 3  Global  Control 

For  closed-loop  control,  we  first  consider  a  discrete-time 
linear  quadratic  regulator.  The  discrete-time  equivalent  of 
the  state  equations,  Eq.  (2-4)  is  (Ref.  11) 

x(/c+  l)-a(*)x(fc)  +  B(fc)u(lt).  fc  =  0,  1.2,... 


(3-7) 


L6 


where  k  denotes  the  time  t~tk  =  kT,  in  which  T  is  the  sampling 
period,  and 


A(fc)-e 


»(*:)- 


aqqt  _  ,  ,  A (k)T  |  Aa(fc)T2  [  ,  [  A \k)T‘ 

1!  2!  i! 

A(k)T2  A2(k)T3  Ai(k)Ti*1 
2!  3!  (t  +  1 ) ! 


B(/c) 


(3-8) 

(3-9) 


are  discrete-time  coefficient  matrices.  Equations  (3-8)  and 
(3-9)  are  obtained  by  regarding  A(0  and  B(<),  Eqs.  (2-5a,b)  , 
as  constant  over  the  steps  kT<t<kT  +  T,  (k  =  0 , 1 , 2 , ...)•  The 
discrete-time  optimal  feedback  control  can  be  written  as 


u(k)  “  %.(k)x(k) ,  k  =  0,1,2,... 


(3-10) 


where  !£(/c)  is  the  discrete-time  control  gain  matrix,  obtained 
by  minimizing  the  discrete-time  quadratic  cost  function  (Ref. 
13) 


N-  1 

J  =  xT(A/)P(A/)x(A/)  +  ]T  [xT  (k)Q(k)x(k)  +  uT  (k)R(k)u(k)]  (3-11) 

fc-0 

The  weighting  matrices  P(A/).Q(0),Q(1) . Q(/V  -  1)  and  R(0),  R(  1 ), 

...,R(jV-1)  are  symmetric.  Moreover,  Q(k)  and  R(A:)are  positive 
definite  and  P(/V)  is  positive  semidefinite,  and  for  our  purposes 
they  are  assumed  to  be  constant  and  diagonal  for  all  N  control 
steps  with  P(A/)  =  Q(A/-  1).  Also,  it  is  assumed  throughout  this 


chapter  that  all  system  states  are  available  for  feedback,  and 
that  there  is  no  estimation  error  involved. 

The  solution  to  this  state-feedback  problem  is  based  on 
the  work  of  either  Kalman  or  Bellman  and  can  be  found  in  a 
variety  of  texts  .  Following  the  procedure  based  on  Bellman's 
Principle  of  Optimality  described  in  Ref.  13,  the  optimal 
feedback  gains  are  calculated  backwards  in  step.  Beginning 
with  i  =  1  and  the  known  P(/V),  we  write 

!*:(N-0  =  -[Sr(A/-0P(N+  1  -i)®(A/-i)+R(A/-Q]"1 

x®T(N-i)P(JV+  l -i)A(N-i)  .  (3-12) 

P(W-i)“[-*(W-0  +  ©(W-0:£(W-0]TP(W+  l  -  OWN-0 

+©(A/-0^(N-0]  +  ^t(N-0R(N-0^(A/-0  +  Q(N-0  (3-13) 

and  continue  until  £(0)  is  calculated.  For  our  system,  the 
final  step  occurs  several  time  steps  after  the  slewing  maneuver 
is  completed,  at  which  time  the  system  can  be  regarded  as 
time-invariant.  Therefore,  the  gains  calculated  soon  after 
the  completion  of  the  maneuver  will  be  treated  as  constant, 
leading  to  steady-state  regulation. 

If  the  system  complexity  results  in  computational 
requirements  in  excess  of  the  controller's  capacity,  given  the 
desired  control  frequency,  then  in  implementing  such  a  control 
law  the  LQR  gains  can  be  calculated  prior  to  maneuver  execution, 
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based  on  the  known  system  parameters  and  slewing  profile,  and 
stored  in  the  computer  for  use  during  and  after  the  maneuver. 
The  disturbance-accommodating  open-loop  control  is  calculated 
during  the  maneuver  and  updated  with  each  time  step. 

In  evaluating  the  utility  of  such  a  control  law  versus 
the  time-invariant  approximation,  we  must  weigh  the  improve¬ 
ments  in  the  time  response  of  the  spacecraft  against  the  added 
computational  complexity  involved  in  a  time-varying  plant. 

3.4  Decentralized  Control 

Implemention  of  a  decentralized  control  law  has  many 
significant  advantages  for  the  type  of  system  under  investi¬ 
gation.  Especially  in  the  case  of  a  spacecraft  with  several 
antennas,  using  many  simple  independent  controllers  reduces 
the  computational  requirements  significantly  compared  to  a 
global  control  system.  The  computations  can  also  be  performed 
completely  "on-line,"  as  opposed  to  a  significant  amount  of 
premaneuver  computations  required  by  global  controls.  This 
computational  savings  could  allow  for  faster,  more  accurate 
sampling  in  the  case  of  discrete-time  control. 

In  this  dissertation,  two  different  types  of  decentralized 
control  are  investigated:  the  first  in  which  individual  points 
on  the  antennas  are  controlled  with  collocated  sensors  and 
actuators  (Collocated  Control)  and  the  second  in  which  indi¬ 
vidual  modes  of  the  independently  modeled  substructures  are 


controlled  (Substructure-Decentralized  Control) .  In  the  case 
of  Collocated  Control,  it  is  assumed  that  both  displacements 
and  velocities  at  the  sensor  locations,  measured  relative  to 
the  appendage  body  axes,  are  available.  In  the  case  of 
Substructure-Decentralized  Control,  each  appendage  has  an 
individual  controller  whose  actuator  forces  are  based  solely 
on  the  generalized  displacement  and  velocity  of  the  substructure 
modes  corresponding  to  that  appendage.  In  both  cases,  the 
rigid-body  motions  of  the  entire  system  are  controlled  inde¬ 
pendently  via  actuators  and  sensors  located  on  the  main 
platform.  As  in  Collocated  Control,  in  the  case  of 
Substructure-Decentralized  Control,  all  states  associated  with 
the  flexible  motion  of  the  antennas  are  available.  Although 
the  controls  are  designed  using  a  decentralized  approach,  the 
response  in  both  cases  is  calculated  by  applying  the  decen¬ 
tralized  control  law  to  the  full  spacecraft. 

3.4.1  Collocated  Control 

3. 4. 1.1  Linear  Control  -  Uniform  Damping 

In  this  case,  a  linear  control  law  is  first  considered 
for  the  controllers  at  all  points.  Control  gains  for  each 
controller  are  designed  and  implemented  using  a  discrete-time 
approach  considering  the  time-varying  nature  of  the  plant. 
Each  of  the  six  controllers  responsible  for  the  rigid  body 
motion  of  the  spacecraft  are  designed  separately,  each  with 
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a  closed-loop  frequency  chosen  well  above  the  fundamental 
frequency  of  the  flexible  appendage  and  damping  equal  to  70% 
of  the  critical  damping.  The  controllers  are  designed  on  the 
basis  of  the  discrete-time  version  of  the  state  equations. 
The  continuous-time  equations  for  the  rigid-body  motions  are 


d_  *<(0 
CLt  X;(0 


Xi(t) 

*,(0 


Bj(Oui(0  > 
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(3-14) 


where 
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0  1 
0  0 


B,(0  = 


0 

ra.'.’CO. 


i  =  1.2 . 6 


(3-  15a, 6) 


in  which  mu  is  mass  or  moment  of  inertia  corresponding  to  the 

rigid  body  mode  i.  Using  full  state  feedback  for  each  controller, 
the  closed-loop  difference  equation  is 

x(fc+ l)-[A(Jt)  +  S(Jt)5C(fc)]x(fc),  Jfe- 0,1,2,...  (3-16) 

where  %,(k)  is  chosen  to  yield  satisfactory  closed-loop  response. 

For  the  point  controllers  located  on  the  appendage(s) , 
the  control  law  is  based  on  the  concept  of  Uniform  Damping 
Control  (Ref.  45) .  This  technique  is  attractive  because  the 
control  design  is  independent  of  the  system  stiffness.  When 
n  discrete  actuators  are  used  on  an  antenna,  the  control  forces 
at  points  Pr  are  given  by 
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Fr(t)  =  -a2MrUr(0-2aMrUr(t),  r-l,2....,/i  (3-17) 

where  U_r(t)  =  u(P r,t)  and  U_r(t)  =  u(P r  ,t)  are  measurements  of 

displacement  and  velocity,  respectively,  at  points  Pr  and 
Mr  ”  I DrP(.F)d  Dr  is  the  mass  associated  with  the  region  on  which 
F_r( t)  acts;  a  is  the  desired  exponential  decay  rate  of  the 
closed-loop  response. 

3.4. 1.2  Nixed  Control 

Next,  we  implement  a  nonlinear  (bang-off-bang)  control 
law  for  the  point  controllers  on  the  appendage.  A  combination 
of  velocity  and  position  feedback  is  used  in  conjunction  with 
a  control  deadzone.  A  dual  elliptical  deadzone  is  used  to 
reduce  control  interaction  between  subsystems  and  slow 
limit-cycle  behavior,  thus  conserving  fuel.  The  two  deadzones 
are  such  that  the  "off"  deadzone  is  smaller  and  lies  completely 
within  the  "on"  deadzone  (Ref.  37)  .  The  "off"  deadzone  is 
defined  by 


<  1 


(3-18) 


(3-19) 
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is  met.  A  state-plane  description  of  the  control  law  is  shown 
in  Fig.  2. 

3. 4. 1.3  Nonlinear  Control 

As  in  the  linear  control  case,  each  rigid-body  mode  is 
considered  separately,  with  its  own  control  law  and  actuator. 
However  in  this  case,  nonlinear  controls  are  used  to  stabilize 
the  rigid-body  motions  as  well.  The  state-plane  switching 
curve  used  for  each  rigid-body  mode  is  that  of  a  time-optimal, 
discrete-time  controller  for  a  double-integrator  system  as 
described  in  Ref.  17.  The  mass  and  inertia  terms,  and  thus 
the  switching  curve  for  each  mode,  are  updated  at  each  control 
step.  However,  because  in  reality  the  system  modes  are  coupled, 
the  resulting  system  control  is  suboptimal.  Figure  3  shows 
these  new  switching  curves  with  control  deadzones  in  the  state 
plane.  The  appendage  controllers  are  of  the  same  type  as  in 
the  previous  case,  in  which  linear  rigid-body  controls  were 
used. 

3.4.2  Substructure-Decentralized  Control 

Next,  we  propose  is  to  control  each  flexible  substructure 
by  means  of  actuator  forces  depending  on  state  measurements 
associated  with  a  given  substructure  alone.  For  example,  in 
the  case  of  a  system  modeled  as  one  or  more  flexible  antennas 
hinged  to  a  rigid  platform,  a  controller  is  implemented  in 
which  each  antenna's  feedback  forces  are  based  on  displacement 
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and  velocity  information  taken  at  discrete  points  on  that 
antenna.  Therefore,  each  antenna  is  modeled  separately  as  a 
cantilever  beam  with  a  substructure  control  law  designed  to 
control  a  finite  number  of  antenna  modes. 

3. 4. 2.1  Linear  Control 

As  in  the  case  of  Collocated  Control,  a  linear  control 
for  both  the  rigid-body  and  the  flexible  motions  is  employed 
first.  The  control  law  used  for  rigid-body  motions  is  unchanged 
from  the  Collocated  Control  cases.  The  equations  of  motion 
for  substructure  e  can  be  expressed  as  the  infinite  set  of 
decoupled  second-order  substructure  modal  equations 

Qi(0  +  a),2Qi(0  =  Qi(0>  t=  1,2,...  (3-20) 


or  as  the  set  of  first-order  state  equations 
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(3-21) 


where  to,  is  the  natural  frequency  and  Q,(0  the  generalized 

control  force,  each  associated  with  the  ith  substructure  mode. 

Of  course,  it  is  neither  necessary  nor  possible  to  control 
all  substructure  modes  using  discrete  controllers.  Hence,  the 
first  step  is  to  truncate  these  equations,  controlling  enough 
modes  to  simulate  the  motion  adequately  and  to  minimize  control 
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spillover  into  the  uncontrolled  substructure  modes.  Using  a 
number  of  actuators  larger  than  or  equal  to  the  number  of 
controlled  modes  helps  reduce  spillover  (Ref.  6) .  Feedback 
gains  are  chosen  for  the  equations  of  the  controlled  sub¬ 
structure  modes  so  that  the  discrete-time  formulation  of  these 
equations  yields  a  damping  factor  and  settling  time  for  the 
fundamental  substructure  mode  of  each  appendage  comparable  to 
the  settling  time  for  the  uniform  damping  case.  Lower  factors 
for  the  higher  modes  are  chosen  so  as  to  yield  a  settling  time 
for  each  equation  comparable  to  that  of  the  fundamental  mode 
for  that  appendage.  The  discrete-time  closed-loop  system 
matrix  can  then  be  constructed  and  the  eigenvalues  checked  at 
key  points  in  the  maneuver  to  assure  convergence.  The  resulting 
r|e  generalized  control  forces  can  then  be  transformed  to  na 
actual  control  forces  for  antenna  e  using  the  transformation 

£.(*)-[*!(-•.)  ♦Kc.j  -  (3-22) 

where  4>g(ze{)  is  the  substructure  participation  matrix,  i.e., 

a  matrix  of  admisssible  functions  evaluated  at  the  location 
of  actuator  e, ,  and  Qe(f) =  [^»i  ^«2 "  'Q«n, jT  is  the  vector  of 
generalized  controls.  If  the  matrix  multiplying  Q„(0  is  not 
square,  then  the  pseudo- inverse,  described  earlier  for  the 
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disturbance  accommodation  control,  should  be  used.  The 
coefficient  matrix  for  the  closed-loop  system  can  then  be 
expressed  as 

•*ei(*)  =  *(*)  +  8(fc)£  (3-23) 


where 


Km  0  Zl  3  0 

0  K22  0  £24 


(3-24) 


in  which  Z\\  and  are  diagonal  6x6  matrices,  the  first  with 

elements  cH(t)  (i  =  1 , 2 , 6)  corresponding  to  discrete-time 
position  feedback  gains  for  the  ith  rigid-body  mode  and  the 

second  with  elements  6ti(0  (t“1.2 . 6)  corresponding  to 

velocity  feedback  gains.  Moreover,  Z22  and  Z2 4  are  the  result 
of  transforming  diagonal  matrices,  with  each  diagonal  element 
corresponding  to  a  feedback  term  from  the  substructure  modal 
eguation  similar  in  form  to  the  rigid-body  equations,  to  the 
configuration  space  using  the  substructure  matrix  of  admissible 
functions;  an  example  is  given  in  detail  in  Appendix  C.  Of 
course,  these  gains  are  constant,  because  they  are  based  on 
the  independent,  time-invariant  substructure,  and  they  are 
zero  for  the  uncontrolled  modes.  The  remaining  blocks  are 
zero  matrices  of  appropriate  dimension. 
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3. 4. 2. 2  Mixed  Control 

The  nonlinear  case  for  Substructure-Decentralized  Control 
is  similar  to  the  Collocated  Control  case  for  controlling  the 
rigid-body  motions.  The  flexible  motion  of  each  antenna  is 
controlled  by  bang-off-bang  controllers  associated  with  sub¬ 
structure  modes.  Because  the  natural  frequencies  of  a  sub¬ 
structure  are  only  indirectly  related  to  the  eigenvalues  of 
the  overall  system,  basing  the  switching  curves  solely  on 
substructure  eigenvalues  would  have  little  value.  Nonlinear 
damping  is  added  to  each  mode  (i.e.,  on-off  control  based  on 
velocity  feedback)  by  means  of  a  modal  force  similar  to  that 
used  in  Ref.  28,  i.e.,  using  the  control  law 

I  -kr.  qr>cr 

/r  =  <  0,  | Qr I  <cr  ,  r  =  1,2,3  (3-25) 

V  Q  r  —  ~  ^  r 

However,  if  the  fundamental  frequency  of  the  system  with 
stabilized  rigid-body  motions  is  so  low  that  settling  time 
requirements  can  be  met  only  through  added  nonlinear  stiffness 
(i.e.,  on-off  control  based  on  position  feedback),  then  a  more 
appropriate  control  law  is 


f-kr,  q  r  ^  d,  u  [qr>-dr  n  qr>cr] 
fr  =  l  0.  t<7r|<dr  n  |qr|  <cr 

V  *r.  qr^-dr  u  [qr<dr  n  qr<-cr] 


r  =  1,2,3  (3-26) 
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In  this  case,  position  feedback  is  used  when  large  displacements 
occur,  while  velocity  feedback  is  employed  for  small  amplitude 
motions.  Figure  4  is  a  state  plane  description  of  this 
criterion. 

3. 4. 2. 3  Nonlinear  control 

The  fully  nonlinear  case  simply  combines  the  nonlinear 
controllers  for  each  appendage  described  in  the  previous  section 
with  the  nonlinear  rigid-body  controller  used  in  the  nonlinear 
collocated  control  case  in  Sec.  3. 4. 1.3. 


4  STOCHASTIC  DISCRETE-TIME  CONTROL 


4 . 1  Introduction 

In  an  effort  to  address  some  of  the  problems  involved  in 
applying  the  previously  discussed  control  techniques,  in  this 
section  I  consider  the  use  of  noisy  actuators  and  sensors. 
In  addition,  I  remove  the  assumption  of  full-state  feedback, 
resulting  in  the  need  to  estimate  unmeasured  states  and  filter 
system  noise  for  the  time-varying  system. 

4.2  Sensors 

Whereas  rigid-body  motions  of  the  main  structure  can  be 
measured  using  conventional  techniques,  such  as  gyros  or  star 
sensors,  measurement  of  elastic  deformations  in  an  antenna 
subject  to  independent  slewing  maneuvers  presents  unique 
difficulties.  This  dissertation  is  concerned  only  with  the 
selection  of  sensors  for  measuring  antenna  elastic  motions. 
Ideally,  these  measurement  devices  should  be  able  to  account 
for  the  changing  nominal  antenna  position  by  either  sensing 
deformations  directly  in  the  antenna  local  coordinate  frame 
or  by  measuring  positions  in  the  rigid-body  frame  and  then 
performing  a  simple  transformation.  Accelerometers,  for 
example,  which  measure  inertial  accelerations,  require  sig¬ 
nificant  processing  of  the  output  to  separate  antenna  local 


28 


29 


accelerations.  This  output,  in  turn,  depends  on  other  measured 
and  estimated  values,  such  as  rigid-body  angular  rotations  and 
rates,  making  accelerometers  a  poor  choice  as  sensors  for 
elastic  motions. 

4.2.1  Strain  Gauges 

Among  sensors  already  used  in  flexible  structure  control 
experiments,  strain  gauges  are  particularly  attractive. 
Indeed,  they  can  measure  antenna  strains,  from  which  dis¬ 
placements  relative  to  the  local  frame  can  be  easily  derived. 
The  strain  gauge  output  voltage  is  proportional  to  the  bending 
strain,  e  of  the  antenna  outer  fibers.  The  beam  strains 
corresponding  to  bending  about  the  x„-  and  y  ,,-axis  are  given 
by 


ex(z.()  =  h5 


<*2Ux(z#(,0 

dzl 


ey(z9() 


32uy(ze(,0 

<?z  I 


1.2 . P„  (4-1) 


where  hx  and  hy  are  the  half-thicknesses  of  the  beam,  whose 

undeformed  neutral  axis  coincides  with  the  local  za  axis,  and 
i  denotes  the  pa  strain  gauge  locations.  Recalling  Eq.  (2-3), 
we  can  write 

<*2ue  (ze,,0 

- j-2 - =  <V(ze()2,(0.  e-  1.2 . N 


(4-2) 
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Assuming  a  circular  cross-section  hx  -  hy  *  h  and  solving  for  the 

vectors  of  elastic  generalized  displacements  of  the  antenna, 
we  obtain 


where 


s.-[*r-(ztl)  *i-(z.,)  - 


(4-4) 


and 


|.*[£*(z.,)  £y(z.,)  £x(Z.2)  «,(Z.  ) 


:(Z*J  £«(Z-..)]r(4-S) 


is  the  vector  of  strain  outputs  from  the  2pa  gauges  located 
on  antenna  e . 

Unfortunately,  strain  gauge  outputs  tend  to  be  noisy,  due 
primarily  to  the  quality  of  the  associated  electronic  equipment 
and  low-output  signal  strength.  Use  of  a  Kalman  filter  to 
compensate  for  this  drawback  will  be  discussed  in  Sec.  4.3. 

4.2.2  Optical  Sensors 

Another  convenient  choice  is  a  photogrammetric  sensor 
(Ref.  49) .  In  this  case,  two  or  more  sensors  are  placed  on 
the  main  spacecraft  for  each  antenna.  Small  infrared  light 


31 


emitting  diodes  (LEDs)  are  then  placed  at  desired  locations 
on  the  antenna  and  fired  sequentially.  As  an  LED  fires,  its 
position  is  sensed  by  a  planar  photodiode.  Combining  data 
from  two  of  these  sensors  and  comparing  this  information  with 
the  known  sensor  and  undeformed  antenna  position  yields  the 
elastic  deformation  at  each  source  point.  One  shortcoming  of 
an  optical  system,  such  as  that  described  above,  is  a  sample-rate 
limitation  that  depends  not  only  on  the  type  of  sensors,  but 
also  on  the  number  of  light  sources.  Because  all  sources  are 
tracked  by  the  same  set  of  sensors,  only  one  source  per 
antenna-sensor  set  can  be  active  at  any  given  instant. 

Another  optical  technique  would  use  the  same  type  of 
sensors,  but  the  optical  signal  would  be  provided  using  a 
combination  of  lasers  mounted  on  the  main  rigid-body  and  small 
mirrors  placed  at  desired  points  on  the  antenna.  Reference 
43  demonstrated  that  such  a  system  can  generate  high  accuracy 
output  at  high  frequencies  (>900Hz) ,  though  its  ability  to 
sense  large  deflections  is  somewhat  limited. 

4.3  Kalman  Filter  Design 

Independent  of  the  measurement  technique  used,  we  assume 
that  the  output  includes  only  elastic  displacements,  and  not 
velocities.  Likewise,  rigid-body  translations  and  rotations 
are  available  while  rates  are  not.  All  measurements  are 
considered  contaminated  by  "white"  noise  whose  statistics  are 
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known.  The  system  itself  is  also  subject  to  white  noise  due 
to  actuator  errors,  unknown  disturbances,  etc.  In  light  of 
this  assumption,  an  estimate  is  required  of  both  the  unmeasured 
and  noisy  states. 

Because  construction  of  separate  substructure  estimators 
is  meaningless,  a  standard  full-order  step-varying  discrete¬ 
time  Kalman  filter  approach  is  proposed  (Refs.  11,13,20).  In 
fact,  unless  a  method  to  directly  measure  elastic  displacements 
and  velocities  is  available,  thus  rendering  additional  pro¬ 
cessing  to  recover  incomplete  or  noisy  states  or  for  compen¬ 
sation  of  the  time-varying  nominal  position  of  the  antennas 
unnecessary,  some  type  of  global  state  estimation  for  feedback 
will  be  required.  Implications  of  a  requirement  for  global 
estimation  on  control  law  design  are  addressed  in  later  sec¬ 
tions. 

Kalman  filter  design  requires  that  the  discrete-time 
equivalent  of  the  system  matrices  A(t)  and  B(l ),  as  well  as  the 
filter  gains,  be  calculated  at  each  time  step.  The  corresponding 
closed-loop  discrete-time  state  equations  are  given  by 

x(k  +  1 )  =  A(k)x(k)  +  S(/c)u(/c)  +  S , (k)w(k) ,  k  =  0,1,2,...  (4-6) 
y(k)  =  Cx(k)  +  v(k)  (4-7) 
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where  k  denotes  the  time  t-tk-kT,  in  which  T  is  the  sampling 
period,  x(k)  is  the  state  vector  and  u(/c)  is  the  control  input 
vector  at  step  k.  Moreover,  w(k )  and  v(k )  are  corresponding 
actuator  and  sensor  noise  vectors,  respectively,  and  A(k)  and 
B(fc)  are  again  the  discrete-time  coefficient  matrices  given 
in  Eqs.  (3-7,8)  .  In  the  case  in  which  rigid-body  displacements 
are  measured  directly  and  strain  gauges  are  used  to  detect 
flexible  motion  in  a  system  with  one  antenna,  C  has  the  form 


^  6X6  0  0  0 

0  5,  0  0 


(4-8) 


where  5,  was  defined  in  Eq.  (4-4). 

The  state  estimate  x(fc)  is  given  by 

x(^)  =  x(fc)  +  ^(t)CT[C^(i)CT^]'1[y(it)-C£(fc)]  (4-9) 

where 

x(k)  =  A(k-  1  )x(k-  1 )  +  S(/c  -  1  )u(/c  -  1)  (4-10) 

is  the  estimate  of  the  state  based  on  the  previous  state 
estimate  and  input, 

M(k)-A(k-  l)2»(Ac-  1  )AT(k  -  1 )  +  S\(k  -  l)QS,(fc-  1)  (4-11) 


and 
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?(k-  l)  =  M(k-  l)-M(k-  l)CT[CM(k-  1  )CT  +  R]~lCM(k-  1  )(4-  12) 


where  M(  0)  is  known  and  R  and  Q  are  the  covariance  matrices 
associated  with  the  sensor  and  actuator  noise  intensities  v(k) 
and  w(k),  respectively. 

For  more  accurate  state  estimation,  the  discrete 
approximation  of  the  known  inertial  disturbance  D(l)cf(0  can 
be  added  to  Eq.  (4-6) ,  so  that 

X(k)  =  A(k-  1  )x(/c -  l)  +  8(k-  1  )u(k-  1  )  +  £>(*:-  l)d(A:-  1)  (4-13) 

where 


D(/c)  = 


A(/c)T2.  A2(/c)73, 
2!  3! 


A‘(A:)T‘° 

(t-1)! 


D(fc) 


(4-14) 


Moreover,  d(Ac)  is  constant  over  each  time  interval  kT  <t  <(k  +  1  )T 

and  is  equal  to  d(!)  evaluated  at  t  =  kT.  Equation  (4-14)  is 
obviated  by  the  use  of  an  open-loop  disturbance-accommodating 
controller,  as  discussed  in  the  next  section. 

4.4  Pseudo-Decentralizad  Control  Law  Design 

As  mentioned  earlier,  when  incomplete  and/or  noisy 
measurements  necessitate  global  estimation,  true  decentralized 
control  is  not  possible.  Instead,  a  pseudo-decentral ized 
control  is  proposed  where  feedback  controls  are  based  on 
estimates  of  antenna  and  rigid-body  motions.  The  system  control 
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laws  are  pseudo-decentralized  in  the  sense  that  substructure 
control  inputs  depend  only  on  state  estimates  for  that  par¬ 
ticular  substructure.  Those  estimates,  however,  are  calculated 
using  a  global  observer  -  in  our  case,  a  Kalman  filter. 

This  dissertation  investigates  two  different  approaches 
to  pseudo-decentralized  control.  As  in  the  case  of  the 
deterministic  system  discussed  earlier,  individual  points 
corresponding  to  actuator  locations  on  the  antennas  are  con¬ 
trolled,  based  on  observer  estimates  of  the  elastic  states  of 
the  points  in  the  antenna  local  coordinate  frame 
(Pseudo-Col located  Control) .  In  the  second  case, 
Substructure-Decentralized  Control  is  employed.  For 
Substructure-Decentralized  Control,  each  appendage  has  an 
individual  controller  whose  actuator  forces  are  based  solely 
on  the  globally  derived  state  estimates  associated  with  that 
particular  appendage.  In  both  cases,  the  rigid-body  motions 
of  the  entire  system  are  again  controlled  independently  via 
actuators  and  sensors  located  on  the  main  platform.  Although 
the  controls  are  designed  by  means  of  a  pseudo-decentralized 
approach,  the  response  is  calculated  by  applying  the  decen¬ 
tralized  control  law  to  the  full  spacecraft. 

Each  of  the  six  controllers  responsible  for  the  rigid-body 
motions  of  the  spacecraft  are  designed  separately,  in  the  same 
fashion  as  in  the  deterministic  cases.  However,  now  the 
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estimated  states  from  the  Kalman  filter  are  used  as  feedback 
for  each  controller.  The  closed-loop  difference  equation  is 
therefore 

x(/c  +  l)  =  ^(fc)x(fc)  +  ®(A:)[5C(A:)r(fc)^|£(/c)],  k  =  0,  1 ,2,...  (4-  15) 

where  u/(fc)  is  the  actuator  noise  vector  and  5C(A:)  is  chosen  to 
yield  satisfactory  closed-loop  response. 

4.4.1  Pseudo-collocated  Control 

The  linear  control  law  for  the  case  in  which  point  * 
controllers  are  located  on  the  appendage (s)  is  again  based  on 
the  concept  of  Uniform  Damping  Control.  For  the  stochastic 
system,  the  uniform  damping  control  forces  at  points  Pe(  are 
given  by 


Fa<ik)--a2Mai0a'(k)-2aMa[0Q'(k),  t-  1,2 . na  (4-16) 


where  f/0((fc)-  u(P0(,k)  and  0_ai(.k)m  ti^Pg^k)  are  estimates  of 


displacement  and  velocity,  respectively,  at  points  P t(  . 
Moreoever,  Q_a  (k)  and  £?#((fc)  are  found  using  the  Kalman  filter 
estimate  of  antenna  states,  qa(k)  and  qa(k).  Hence, 


0a(k)m<t>a(Pai)qa(k) 


e  =  1 , 2 1 V  (4-17) 


3: 


On  the  other  hand,  Mi(  and  a  retain  the  same  meaning  as  in  the 
deterministic  case. 

Next,  a  nonlinear,  bang-off-bang  control  law  is  implemented 
for  the  point  controllers  on  the  appendage  using  the  same  type 
of  state-plane  switching  curve  as  in  Section  3.4.1.  Because 
the  controls  are  now  based  on  estimates  of  the  local  position 
and  velocity,  the  deadzone  "on"  and  "off"  criteria  become 


(4-18) 


(4-19) 


For  the  fully  nonlinear  case,  the  rigid-body  control  is 
the  same  as  that  used  in  the  deterministic  case  except,  of 
course,  that  now  the  nonlinear  controls  are  based  on  the 
estimated  rigid-body  states. 

4.4.2  Substructure-Decentralized  Control 

Next,  the  system  is  controlled  in  a  manner  similar  to  the 
earlier  Substructure-Decentralized  Control  (Sec.  3.4.2), 
except  that  the  feedback  control  for  a  given  flexible  sub¬ 
structure  is  based  on  estimated  states  associated  with  that 
substructure  as  given  by  a  global  Kalman  filter.  Otherwise, 
the  same  control  laws  can  be  used  for  the  linear,  mixed  and 
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nonlinear  cases.  For  convenience,  the  term  Substructure- 
Decentralized  Control  is  retained  although  the  technique  is 
actually  a  pseudo-decentral ized  control. 


5  NUMERICAL  EXAMPLES 


5 . 1  Introduction 

The  mathematical  model  considered  consists  of  a  rigid 
platform  with  one  flexible  antenna,  represented  by  a  uniform 
beam,  hinged  to  the  platform  at  one  end  and  free  at  the  other 
(Fig.  5)  .  The  antenna  is  slewed  through  a  45°  angle  about  the 
x ,,-axis.  The  beam  is  discretized  in  space  using  five  admissible 
functions  for  each  displacement  component.  The  admissible 
functions  for  the  y  ..-displacement  component  are  taken  in  the 
form  of  cantilever  modes 

<()y  =  -(cos(3yz-cosh  (3  y  z )  +  C/(sin(3/z-sinh|3/z),  j  -  1,2, 3, 4, 5 

This  is  consistent  with  the  fact  that  the  flexible  antenna 
represents  a  slewing  cantilever  beam,  where  the  slewing  angle, 
is  a  given  function  of  time.  The  admissible  functions  for  the 
x ,-displacement  component  have  the  same  form.  The  platform 
has  one  actuator  corresponding  to  each  rigid-body  degree-of- 
freedom  and  the  antenna  has  three  sets  of  bi-directional 
actuators  located  at  z,  -  lt/3,2lt/3,  l,.  The  following  data  was 
used  in  the  computer  simulation: 
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Structural  damping  equal  to  0.2%  of  critical  was  added  to  the 
simulation  model  for  each  of  the  system  elastic  modes  in  the 
premaneuver  configuration. 

The  time-optimal  slewing  profile  is  shown  in  Fig.  6.  For 
all  cases,  the  angular  acceleration  rate  of  the  antenna  was 
±n/8  rad/s,  so  that  the  entire  maneuver  is  completed  in  about 
2.8  s.  The  uncontrolled  motion  of  the  platform  relative  to 
the  inertial  space  and  the  elastic  displacement  of  the  antenna 
tip  are  shown  in  Fig.  7. 

When  only  the  rigid-body  modes  are  controlled  linearly, 
the  resulting  deterministic  uncontrolled  tip  displacement  of 
the  antenna  in  the  y#-direction  is  quite  large,  as  can  be  seen 
in  Fig.  8. 

5.2  Deterministic  Control 

In  the  case  of  control  using  the  step-varying  LQR,  the 
weighting  matrices  P  =  Q  =  100  x  [| ]24x24  and  R  =  .001  x  [|]12xl2  were 
used.  Figure  9  shows  the  tip  displacement  for  the  case  in 
which  open-loop  disturbance  accommodation  is  employed  together 


with  the  LQR.  Figure  10  shows  the  same  information  without 
disturbance  accommodation.  In  all  cases,  except  for  full 
nonlinear  control,  the  rigid-body  displacements  were  insig¬ 
nificant  when  compared  to  the  elastic  displacements  of  the 
antenna . 

Linear  Collocated  Control  using  Uniform  Damping  was 
implemented  using  a  decay  rate  a =5.0.  The  rigid-body  motion 
was  controlled  using  a  linear  control  law  with  uo„  =  10 rad/s  and 
£  =  0.70.  Figures  11  and  12  show  displacements  with  and  without 
the  open-loop  disturbance  accommodating  control,  respectively. 

For  the  case  of  nonlinear  antenna  control  using  the 
parameters 

/ ,  =  0. 18A/  f  2  -  0.5 A/  f  3  =  0.22 /V 
cr=.05m/s,  aoffr  =  -00 1  m/s  aon=.01m/s 
b  off=  .003m/ s  bon=.03m/s,  r=  1,2,3 

the  resulting  tip  displacement  is  as  shown  in  Figure  13.  Figures 
14-16  display  the  time  history  of  the  substructure  actuator 
forces. 

In  the  case  of  nonlinear  control  of  both  the  rigid-body 
and  the  elastic  motions,  the  parameters  used  were  as  follows: 
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/ !  =  0. 18N  f  z  =  0.5/V  /3  =  0.22N 

cr  =  .04m/s  aoff=.0\m/s  aon=.02m/s 
boff=.05m/s  bon  =  AOm/s,  r=  1,2,3 
Fx  =  Fy  =  Fz  =  0.8/V  Mx  =  My  =  Mz  =  1.6N 
a3i  =  a^  =  az  =  .001  m/s  a6^~  a6^  =  =  .001  rad/  s 

b  x  =  b  y  =  b  z=  .003m  bQ^  =  b9^  =  60j  =  .003 rad 

The  corresponding  antenna  displacement  and  actuator  histories 
are  shown  in  Figs.  17-20,  Because  the  most  significant 
rigid-body  motion  is  rotation  in  the  plane  of  the  maneuver, 
Figs.  21  and  22  show  the  in-plane  rotation  and  the  actuator 
torque,  respectively. 

In  the  case  of  linear  Substructure-Decentralized  Control, 
the  discrete-time  control  gains  for  each  of  the  three  controlled 
substructure  modes  were  designed  so  as  to  yield  an  exponential 
decay  rate  of  a  =  5.0.  The  antenna  tip  displacement  with  and 
without  disturbance  accommodation  is  shown  in  Figs.  23  and  24, 
respectively. 

A  nonlinear  control  next  replaces  the  linear 
Substructure-Decentralized  Control  with  the  following 
switching  parameters: 

c,  =  c2  =  c3  =  .005m/s ,  k ,  =  0.80 N  ,  k2  =  k3  =  0.40 N 
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Nonlinear  damping,  as  described  earlier,  was  used  for  the  first 
substructure  mode  with  di~  0.01m.  The  resulting  tip  dis¬ 
placements  are  shown  in  Figure  25  and  the  actuator  force 
histories  are  shown  in  Figures  26-28. 

In  the  case  in  which  both  the  rigid-body  and  the  elastic 
control  is  nonlinear,  use  of  the  parameters 

k ,  =  0.80N  k2  =  k3  =  0.40N 

c1  =  c2  =  c3  =  .005m /s  d ,  =  0.0 1  m 

Fx  =  Fy  =  Fz  =  \.0N  Mx  =  My  =  MZ  =  2.0N 

ax  =  ay  =  az-  .001  m/s  a6x  =  a6^  =  a6^  -  .001  rad/ s 

b  x-  b  y  =  b  2  =  .003m  b6*  =  b0^  =  60j  =  .001  rad 

resulted  in  the  antenna  tip  displacements  and  actuator  histories 
shown  in  Figs.  29-32.  Figures  33  and  34  again  display  the 
rigid-body  in-plane  rotation  and  actuator  torque. 

5.3  Stochastic  Control 

The  Kalman  filter  was  based  on  a  truncated  model  of  the 
simulated  system.  In  this  case,  the  antenna  was  modeled  by 
only  three  admissible  functions  in  each  direction,  corre¬ 
sponding  to  the  three  controlled  substructure  modes.  It  was 
also  assumed  that  the  actuator  noise  manifests  itself  as  unknown 
forces  and  torques  on  the  system.  Therefore,  we  choose 
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L 1  12x12  J 

In  the  simulation,  "white”  noise  was  added  to  all  outputs, 
y(/c),  at  each  sampling  time  to  simulate  sensor  noise.  A  similar, 
uncorrelated  "white”  noise  was  added  to  each  system  state 
(including  uncontrolled  flexible  states)  at  every  simulation 
time  step  to  simulate  actuator  noise  and  unknown  disturbances. 
The  "noise  level"  used  for  each  term  was 
£{ot2(/c)>  »  (fc)>  =  10'4.  For  simplicity,  the  sensor  noise 

covariance  matrix  R  is  taken  as  an  identity  matrix  of  appropriate 
dimension  and 
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Linear  Pseudo-Collocated  Control  using  uniform  damping 
was  implemented  first,  again  using  a  decay  rate  a =5.0.  The 
rigid-body  motion  was  also  controlled  using  linear  control 
with  (A)  =10  rad/s  and  £  =  0.70,  i.e.,  the  same  as  in  the 
deterministic  case.  Figure  35  shows  the  antenna  tip  dis¬ 
placement  with  and  without  the  open-loop  disturbance  accom¬ 
modating  control,  respectively. 

Using  the  same  parameters  ki,cl,al,b c,  the  resulting  tip 
displacement  is  as  shown  in  Fig.  36.  For  comparison,  the  case 
with  only  rigid-body  control  (Fig.  8)  is  also  shown.  Figures 
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37-39  display  the  time  history  of  the  substructure  actuator 
forces. 

For  nonlinear  control  of  both  the  rigid-body  and  the 
elastic  motions,  the  parameters  used  were  as  follows: 

/j=0.18iV  /2  =  0.5A/  f3  =  0.22N 

cr  =  .04m/ s  a0/f  =  .Olm/s  aon=.02m/s 
boff=.03m/s  bon=AOm/s,  r=  1,2,3 
F  x  =  F  y  =  F  z  =  0.8N  Mx=My  =  Mz  =  1.6  A/ 
a*  =  ay  =  a. z=  .005 m/s  =  a6*  =  a6^  =  . 005  rad /s 

bx  =  by  =  bz=  .015m  bex  =  bey  =  be ,  =  -015  rad 

The  resulting  antenna  tip  displacements,  antenna  actuator 
forces,  in-plane  rotation  and  actuator  torque  are  shown  in 
Figs.  40-45,  respectively. 

In  the  case  of  linear  Substructure-Decentralized  Control, 
the  discrete-time  control  gains  for  each  of  the  three  controlled 
substructure  modes  were  also  chosen  identical  to  those  in  the 
deterministic  case.  The  antenna  tip  displacement  with  and 
without  disturbance  accommodation  is  shown  in  Fig.  46.  In 
Figs.  46-56,  the  term  Substructure-Decentralized  Control 
(Pseudo)  is  used  to  distinguish  between  these  and  those  from 
the  deterministic  cases. 


The  linear  Substructure-Decentralized  Control  was  then 
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replaced  with  the  nonlinear  control  used  in  Sec.  5.2.  The 
resulting  tip  displacement  is  shown  in  Fig.  47,  again  compared 
with  rigid-body  control  alone  (Fig.  8) ;  the  actual  actuator 
force  histories  are  shown  in  Figs.  48-50. 

For  the  case  in  which  both  the  rigid-body  and  the  elastic 
control  is  nonlinear,  use  of  the  same  parameters  as  in  the 
deterministic  case  yields  the  results  displayed  in  Figs.  51-56. 


6  CONCLUSIONS 


This  dissertation  considers  the  control  of  a  spacecraft 
with  retargetable  flexible  antennas  using  discrete-time 
techniques.  It  significantly  expands  upon  earlier  work  that 
derived  the  continuous-time  equations  of  motion  for  this  type 
of  systems  and  applied  global  control  laws  based  on  a  noise-free, 
full  state  feedback  system  that  was  considered  largely  time- 
invariant.  Chapter  3  of  this  dissertation  presents  a 
discrete-time  approach  for  designing  an  adequate  control  law 
that  accounts  for  the  time-varying  nature  of  the  system,  but 
still  assumes  full  state  feedback  and  no  sensor  or  actuator 
noise.  First,  a  global  approach  using  a  discrete-time, 
step-varying  linear  quadratic  regulator  was  implemented  to 
control  the  time-varying  mathematical  model.  Because  a  global 
controller  requires  extensive  computations,  the  step-varying 
feedback  gains  can  be  precalculated  and  stored  for  use  during 
a  known  spacecraft  maneuver  of  a  system  with  many  appendages. 
Next,  two  decentralized  approaches  were  proposed:  Collocated 
Control  and  a  new  method  designated  Substructure-Decentralized 
Control.  Both  offer  the  advantage  of  lower  computational 
effort,  while  still  accounting  for  the  time-varying  nature  of 
the  system.  Whereas  Collocated  Control  suppresses  the  elastic 
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motion  at  specific  points,  Substructure-Decentralized  Control 
suppresses  the  elastic  appendage  substructure  modes,  taken  as 
cantilever  modes  in  the  case  of  a  beam-like  antenna.  In  each 
of  the  decentralized  cases,  linear  and  nonlinear  (on-off) 
schemes  were  presented  for  controlling  the  elastic  motion  of 
the  antennas  and  the  rigid-body  motions  of  the  system.  The 
linear  Collocated  Control  law  was  based  on  a  discrete-time 
application  of  Uniform  Damping.  Note  that  on-off  controls 
used  to  suppress  elastic  motions  translate  into  quantized 
actual  controls  in  the  Substructure-Decentralized  Control  case. 
The  nonlinear  control  of  the  rigid-body  motions  assumed  full 
decoupling  of  the  rigid-body  modes  and  used  a  switching  curve 
based  on  a  discrete-time,  step-varying  interpretation  of  a 
time-optimal  controller  for  a  double-integrator  system.  In 
all  global  (LQR)  and  decentralized  cases  using  fully  linear 
control,  a  discrete-time  disturbance-accommodating,  control 
designed  to  counteract  the  inertial  disturbances  caused  by  the 
known  antenna  maneuver  profile  was  also  included. 

Chapter  4  considers  the  effects  of  noisy  actuators  and 
sensors,  as  well  as  the  reality  of  incomplete  state  feedback, 
was  addressed.  A  step-varying,  full-order  discrete-time  Kalman 
filter  was  developed  to  estimate  the  system  states  upon  for 
feedback  control.  Strain  gauges  were  selected  to  measure 
•elastic  deformation  in  the  slewing  antennas  because  of  the 
small  amount  of  processing  required  to  determine  system  states 


from  their  output.  Pseudo-decentralized  control  was  then 
proposed  whereby  control  forces  used  to  suppress  the  elastic 
motion  of  an  individual  antenna  are  based  on  the  globally 
derived  estimates  of  states  corresponding  to  that  specific 
antenna.  Indeed,  if  any  processing  of  antenna  sensor  outputs 
is  required  to  account  for  the  relative  motion  of  the  antenna 
local  coordinate  frame  during  a  maneuver,  true  direct  output 
feedback  control  with  collocated  sensors  and  actuators  is  not 
possible.  Similarly,  a  requirement  of  state  estimation  due 
to  noise  of  incomplete  measurements  would  make  true  decen¬ 
tralized  control  in  general  impossible. 

The  decentralized  controllers  applied  to  the  deterministic 
system  were  therefore  selected  for  use  on  the  stochastic  system. 
The  difference  here  is  that  control  forces  are  based  on  estimates 
of  substructure  modal  states  in  the  case  of  Substructure- 
Decentralized  Control  and  estimates  of  the  position  and  velocity 
of  specific  points  on  the  antenna,  measured  in  the  antenna 
local  coordinate  frame,  in  the  case  of  pseudo-col located 
control . 

A  numerical  example  using  a  deterministic  system  with  one 
antenna  was  presented  and  the  results  of  each  case  compared. 
In  each  of  the  linear  cases,  the  spacecraft  was  controlled 
quite  well,  with  both  decentralized  methods  offering  per¬ 
formance  comparable  to  the  more  computationally  intensive 
step-varying  LQR.  While  the  linear  controls  performed 
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significantly  better  than  the  mixed  or  nonlinear  cases,  as 
might  be  expected,  the  mixed  controller  methods  performed 
adequately,  reducing  tip  motion  for  the  5  m  antenna  to  under 
4  mm  within  one  second  of  maneuver  completion,  with  minimal 
limit-cycle  behavior.  The  fully  nonlinear  controller  case 
using  Substructure-Decentralized  Control  also  performed  well, 
converging  quickly  with  little  noticeable  limit-cycle  behavior. 
The  Collocated  Control  case  exhibited  more  limit-cycle 
behavior. 

Next,  numerical  examples  using  a  stochastic  version  of 
the  same  system  with  one  antenna  was  presented.  Strain  gauges 
to  sense  elastic  motion  and  a  step-varying  Kalman  filter  were 
used  and  the  results  for  each  different  control  law  again 
displayed.  Both  linear  pseudo-decentral ized  methods  controlled 
the  spacecraft  quite  well,  with  performance  comparable  to 
earlier  results  based  on  noise-free  full-state  feedback.  Once 
again,  the  linear  controls  performed  much  better  than  the  mixed 
and  nonlinear  ones.  Both  mixed  control  methods  performed 
adequately,  comparable  to  the  determistic  applications  of  the 
same  controllers.  The  fully  nonlinear  controllers  performance 
was  also  similar  to  that  of  their  deterministic  counterparts. 

In  conclusion,  simulations  confirm  that  for  both  the 
deterministic  and  stochastic  case,  linear  and  nonlinear 
discrete-time  antenna  controls  based  on  Substructure- 
Decentralized  Control  are  able  to  suppress  elastic  motions  of 
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a  slewing  appendage  in  a  time-varying  spacecraft  effectively. 
Strain  gauges  can  conveniently  sense  elastic  states  associated 
with  each  appendage  and,  when  coupled  with  a  global  step-varying 
estimator  (or  Kalman  filter  for  a  noisy  system) ,  can  provide 
feedback  for  pseudo-decentralized  substructure-decentralized 
controllers.  Based  on  globally  estimated  states,  each  con¬ 
troller,  in  turn,  effectively  controls  the  elastic  motions  of 
its  individual  antenna. 
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Fig.  1  Rigid  Platform  with  Flexible  Appendages 


Fig.  5  Mathematical  Model 


Angular  Displacement  (rad) 
Angular  Velocity  (rad/s) 


63 


Displacement  (m;rad) 


6 


Displacement  (m 


66 


Fig.  9  Tip  Displacement  for  LQR  with  Disturbance 

Accommodation 
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Fig.  11  Tip  Displacement  for  Linear  Collocated 
Control  with  Disturbance  Accommodation 
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Fig.  14 
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Fig.  15 
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Fig.  16  Antenna  Actuator  Force  for  Linear  Rigid-Body 
Control,  Nonlinear  Collocated  Antenna  Control,  z 3=£a 
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Fig.  17  Tip  Displacement  for  Nonlinear  Rigid-Body  and 
Collocated  Antenna  Control 
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Fig.  23  Tip  Displacement  for  Linear  Substructure- 
Decentralized  Control  with  Disturbance  Accommodation 
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Fig.  24  Tip  Displacement  for  Linear  Substructure- 
Decentralized  Control  without  Disturbance  Acci.  jnodation 


0.00 


2.00 


4.00 


6.00 


Time  (s) 


Fig.  25  Tip  Displacement  for  Linear  Rigid-Body 
Control,  Nonlinear  Substructure-Decentralized  Control 
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Fig.  29  Tip  Displacement  for  Nonlinear  Rigid-Body 
and  Substructure-Decentralized  Antenna  Control 
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Fig.  35  Tip  Displacement  for  Linear  Pseudo-Collocated 
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Fig.  36  Tip  Displacement  for  Linear  Rigid-Body  Control, 
Nonlinear  Pseudo-Collocated  Antenna  Control 
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Fig.  40  Tip  Displacement  for  Nonlinear  Rigid-Body  and 
Pseudo-Collocated  Antenna  Control 
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Fig.  47  Tip  Displacement  for  Linear  Rigid-Body  Control 
Nonlinear  Substructure-Decentralized  Antenna  Control 

(Pseudo) 
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Fig.  56  x-Axis  Actuator  Torque  for  Nonlinear  Rigid-Body 
and  Substructure-Decentralized  Antenna  Control  (Pseudo) 


APPENDIX  A 


If  a(0  is  a  three-dimensional  vector,  i.e., 

a(.t)  =  [al(t)a2(t)a3(t)]T ,  then  the  skew  symmetric  matrix  <5(0/ 
called  the  vectrice  (Ref.  14)  is  defined  as 


<5(0  = 


0 

a3(0 

-a2(0 


~a3(0 

0 

a,(0 


a2(0 

-a,(0 

0 


(/l-l) 


Using  this  expression,  the  cross  product  of  two  vectors  is 
then  given  by 


a(0><r  =  <3(0r  =  -fa(0  ( -  2) 

Furthermore,  the  skew  symmetric  matrix  of  the  cross-product 
is 


[<56]  =  (56-6<5 


M-3) 


and  the  time  derivative  of  a  direction  cosine  matrix,  for 
example 


C(0  = 


C02C03  C0,S03+S01S02C03 

- C02S03  C0|C03  —  S0jS02S03 

s02  - S0 , C02 


S0,S03-C0,S02C03 
S0 , C03  + C0J  S02S03 
C0, C02 


(A-4) 


where  c0t*cos0it  s0(  =  sin0„  can  be  expressed  as 
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C  =  -(JoC  (AS) 

or 

CT  •*  CTco  (A- 6) 

where 

O)  =  ©0  ( /!  -  7) 

and 

c02c03  s03  0~ 

©(0=  -c02s03  c03  0  ( /l  -  8) 

s02  0  1 


APPENDIX  B  MATRIX  EXPRESSIONS 


The  components  of  the  system  state  equations  (Eq.  2 -5a) , 
M(t).  K(i)  and  G(0  have  the  form 
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where 


n  *  s 

m(  =  mr+£ms  Sf-(  r(dm(  5,  =  V  (m8/*0#  +  Ej5eEe) 

«* 1  J  m,  a- 1 
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(5-4, 5,6) 

(5-7,8) 

(5-9) 


M,=  ]  4>a=  <t>sd m#  $#=  j  f,<t>edm,  (5-  10,  1 1 , 12) 

i/m,  ./  m, 

fi.(^4)  =  f  <*>T*6>a<i>'<irna  H,(u >  )-  f  1>[(jo^,dm(  (5-13,14) 

f  *<e  =  [<l>e,<i>e]  (5-15,16) 


in  which  the  symbol  [  ,  ]  represents  an  energy  inner  product 
(Ref.  24) . 


G22=  £(El(2^ala-trlscus)Ee  +  270tEl[5eu)e]Ee)  (5-17) 

«-  1 

G23  =  270#E[coe5e+E[a)e<fe  +  E[je(u)e)  (5-18) 

The  matrix  B'(<)  (Eq.  2-5b)  is  given  by 
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The  inertial  disturbance  vector  is  defined  as 
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APPENDIX  C  8 DC  LINEAR  FEEDBACK  MATRICES 


The  submatrices  (from  Eq.  3-23)  associated  with  the 
feedback  gains  for  linear  Substructure-Decentralized  Control 
applied  to  the  flexible  appendages  of  the  system  in  question 
are  given  by 

«»c*)-[+r(u) 

where,  for  example,  in  the  case  of  one  antenna. 
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Note  that  Cyu  =  byu  =  0  if  the  ith  substructure  mode  in  the 

y .-direction  is  uncontrolled,  and  the  same  can  be  said  for 
substructure  modes  in  the  x, -direction. 
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